Exactly soluble spin-i models on three-dimensional lattices and non-abelian statistics 

of closed string excitations 
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Exactly soluble spin-| models on three-dimensional lattices are proposed by generalizing Kitaev 
model on honeycomb lattice to three dimensions with proper periodic boundary conditions. The 
simplest example is spins on a diamond lattice which is exactly soluble. The ground state sector 
of the model may be mapped into a p-wave paired state on cubic lattice. We observe for the first 
time a topological phase transition from a gapless phase to a gapped phase in an exactly soluble 
spin model. Furthermore, the gapless phase can not be gapped by a perturbation breaking the time 
reversal symmetry. Unknotted and unlinked Wilson loops arise as eigen excitations, which may 
evolute into linked and knotted loop excitations. We show that these closed string excitations obey 
abelian statistics in the gapped phase and non-abelian statistics in the gapless phase. 

PACS numbers: 75.10. Jm,03.67.Pp,71.10.Pm 



Introductions: Non-abelian anyons in two dimensions 
provide a promising candidate for quantum computation, 
which is topologically protected from decoherence 
Recently, this topic has attracted great interests 0] ■ The 
most possible system in which the quasiparticles obey 
non-abelian statistics is two-dimensional electrongas in 
a fractional quantum Hall state with v = b/2 [1, |J. 

In search for non-abelian statistics in exact soluble 
models, Kitaev proposed an exactly soluble spin model 
on honeycomb lattice and showed that in a wider pa- 
rameter region, the vortex excitations obey non-abelian 
statistics [5]. Since the exotic statistics of the vortices, 
Kitaev model has attracted many research interests 0| ■ 
It has been shown that Kitaev model is in the same uni- 
versality class with px + ipy-wa.ve paired state Q . 

The point-like anyons are strictly restricted in two di- 
mensions. In three dimensions (3-d), the point-like par- 
ticles can only be either bosons or fermions. However, 
the exotic statistics may arise in the closed string exci- 
tations. In a seminal paper 9], Witten has shown that 
the Wilson loops in 3-d Chern-Simons field theory obey 
non-abelian statistics, which is closely related to the mon- 
odromy matrix in conformal field theory [ll| . The statis- 
tics of unknotted, unlinked closed strings is described by 
loop braid groups The quantum loop gas applied to 
topological quantum computation is a rapidly developing 
field 13,] . Topological quantum order in brane systems 
was also an interesting subject [l4| . 

In this paper, we will generalize Kitaev model to that 
on a diamond lattice as well as multi-layer honeycomb 
lattices in which the closed string excitations obey non- 
abelian statistics. This 3-d generalization of Kitaev 
model is also exactly soluble and the ground state sector 
is equivalent to 3-d p-wave paired state. The phase di- 
agram consists of two topological phases, a gapped one 
and gapless one. The gapped phase is a strong pairing 
phase whose topological nature is characterized by Hopf 
invariant and then the ground state is topologically triv- 
ial. The singularity of the Hopf mapping in the gapless 
phase implies that the ground state is topologically non- 



FIG. 1: The diamond and cubic lattices, (a) The diamond 
lattice and four types of links, (i=l,2,3) corresponding to 
X-, y-, and z-links and ez to ti-links. (b) The equivalent cubic 
lattice. 

trivial. Removing the singular points, a non-zero wind- 
ing number may appear. Thus, in the gapped phase, 
the closed strings obey abelian anyonic statistics while in 
the gapless phase, the strings obey non-abelian statistics, 
which is tantamount to Wilson loops in SU(2)2 Chern- 
Simons field theory Q. The Majorana fermion excita- 
tions in the gapless phase is always gapless even there is 
a perturbation with time-reversal symmetry breaking. A 
topological phase transition between gapped and gapless 
phases [1(1] is found for the first time in an exactly soluble 
model. 

Model and Solution: We begin with a Kitaev-type cou- 
pled spin-i model on a diamond lattice (Fig. [IJi), whose 
Hamiltonian reads 



i+e2 



T ^^^^ 



(1) 



where the translational invariant vector is i = mei -|- 
ne2 + le^ — (m + n + Vje^- This Hamiltonian has a set 
of integrals of motion {Sp,p = 1,2,...} in which anyone 
commutes with the Hamiltonian as well as another mem- 
ber in the set. Sp is a loop operator, i.e., an unknotted 
and unlinked closed string operator along e^3 = 62 — ei if 
we take periodic boundary condition along ej^3. The ex- 
plicit form of Sp is Sp = 

reduces to S'„ ^ 
strain ai 



^i+ne 



me±3 "i+ei'^i i "i+e2 

nm=o ^"^me^a Under the periodic con- 
. Each Sp projects a sector of the to- 
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tal eigenspace of the Hamntonian Hd since [Sp, Hd] = 0. 
This guarantees the existence of exact solutions. 

Expressing spins by Majorana fermions [H], one has 
= ibxC, (T^ = ibyC, and = ibzC with constraint 
bxbybzC — 1. Defining u 



ibi '^b^ for the x- and y- 



links and Uij — ib^bj for z- and w-links, the Hamiltonian 



Hd reads 



Hd = 



s 
s 

-}- iJx ^ ^ '^s,bw{ps.b^s—e\.w *^s,i(;Cs — ei ,6) 
s 

+ Jy Sz Jy partners, (2) 



where s is the position of the z-hnks. Note that 
[Hd.Uij] = 0, ufj — 1 and Sp is the product of a set 
of Uij. Therefore, one may take Uij = ±1 and the eigen 
value of Sp may be determined by product of Uij . Ac- 
cording to Lieb's theorem [l^, the ground state is in- 
cluded in the sector with all Ug^bw = 1, which is con- 
sistent with all Sp = 1. If there are odd numbers of 
Ubw = ^1 along a loop, this gives Sp = —1. This 
is a loop excitation and in fact is a nontrivial Wil- 
son loop. Defining the link fermions living in 2-links, 
ds = (csM + jcs,,u)/2, dl = {cs,b - ics,w)/'2. and deform- 
ing the lattice to a cubic lattice (Fig. ^p) in a similar 
way deforming the honeycomb to a square lattice in two 
dimensions [J Q , the loop excitation-free Hamiltonian 
may be written as 

Ai,, 



+ i- 



(3) 



where dp is the Fourier component of ds ; the dispersion 
and the pairing functions are 

= Jz - Jx COSPx - Jy COS Py - Jy cospz, (4) 

Aa,p = Aax sinpx + Aay smpy + Aaz sinpz, a =1,2. 

This is a p-wave paired state. At present, Aib = 
and A26 = Jb with b = x,y,v. Furthermore, 
one can follow a previous work by Wang and one of 
the authors (Y.Y) [8] and introduce three-spin cou- 
plings and four-spin couplings etc. This leads to the 
parameters Jx^y and Aa-x,y become tunable. To let 
Jz and Aaz be tunable, one needs to add terms like 

<b<-He,,b + <«,c^J+e,,™ and <,bcrj+e,,„- Thesc are 
not nearest neighbor couplings. One can rewritten them 

as '^s,biK+e^,wK+e^,w)'^s+e^,b^ '^twiK,bK,b)'^s+e^.w ^ud 
'^s,w{K,bK,b)iK+e,,wK+e,,wys+e,,^,b bcCaUSC bzbz = 1. 

Thus, these terms can also be bilinear since My are cor- 
rectly inserted. Thus, the model with these terms added 
is still exactly soluble but these pairing parameters Aab 



become tunable. Hereafter, we discuss this general p- 
wave paired state. The quasiparticle excitations are gov- 
erned by the BdG equations 



EpUp = ^pUp - A* Dp, EpVp = ~^pV 



ApUp 



(5) 



where Ep = 



p , (Ai_p)2 + (A2,p)2 is the quasiparti- 
cle dispersion, Ap = Ai p+iA2,p, and {up, Vp) are the co- 
herence factors with |upp = 5(l + ^),|wp|^ = 5(1 ~"E^) 
and Vp/up = -{Ep - Cp)/A;. 

The phase diagram can be drawn in a similar way to 
that for the two-dimensional copy [3] . The phase bound- 
ary is determined by ^p = for any Aab- That is, the 
band insulator/free Fermi gas transition determines the 
phase boundary, which is given by |cosp*| = 1 with 
p* = (0, 0, 0), (0, 0, ±7r),(0, 0),(±7r, 0, 0), (0, ±7r, ±7r), 
(±7r, 0, ±7r),(±7r, ±7r, 0) and (±7r, ±7r, ±7r) where Jz ± 
Jx ± Jy ± Jv — 0. The phase where ^p > is always 
gapped. In the p-wave sense, this is the strong pairing 
phase as that in two dimensions [l^. However, out of the 
gapped phase, there are a pair of p* so that i?p. = for 
general A^fc. This means that there is a phase transition 
from the gapped phase to a gapless phase even the time 
reversal symmetry (or the generalized inversion symme- 
try [s^l) is broken. 

Topology of Phases: Now, the question is that whether 
the phase transition from the gapped to gapless phases 
is a topological phase transition. First, we do not ob- 
serve the spontaneous breaking of any continuous sym- 
metry. This implies the phase transition may be topolog- 
ical. However, in general, a topological phase transition 
requires an energy gap between the degenerate ground 
state and any excitation state. To check if the phase 
with a gapless Majorana fermion excitation is topologi- 
cally non-trivial, we calculate the topological invariants 
of both phases in the continuum limit. The momentum 
space is x S^ due to the periodic boundary condition. 
Note that the coherent function ip'' = {u*,v*) defines a 
mapping from p G D"^ x S^ to {u, v) £ 5^. The unit vec- 
tor m = = (ReAp, — ImAp, ^p)/_Ep parameterizes 
this mapping. If we compact x S^ to S^ by defin- 
ing m = mo on the boundary torus T^, the mapping is 
called the Hopf mapping. Associated with Hopf map- 
ping, there is a topological invariant, the Hopf invariant, 
which is defined by 



(6) 



with Ai{p) = |(-0t5p.V - i'dp^i'^) and Fij{p) = 
dp-Aj{p) — dpjAi{p) = m ■ {dp.m x dp^m). This is an 
abelian Chern-Simons in the momentum space. On the 
other hand, we know that the Hopf invariant describes 
the linking numbers of closed strings Due to the 

duality between the cubic lattice and dual momentum 
lattice, the linking number in momentum space is equal 
to that in the co-ordinate space. For the strong pairing 
ground state, this linking number is zero because there is 





(a) (b) 



FIG. 2: (a) The bonding configuration of the multilayer model 
Hz- (b) The bonding configuration of another multilayer 
model mentioned in [23j. 

no closed string excitation. One may also directly prove 
that {u{oo),v{(xi)) = (m(0), ^(O)) = (0,1) in the ground 
state and thus, Hf=0. 

Note that the unit vector m(p*) is singular in the gap- 
less phase because at Dirac points p*, Ep* = 0. For a set 
of given parameters in the gapless phase, there are two 
such Dirac points except at the phase boundary where 
Pi = P2 ~ 0. The unit vector m can not be defined in 
these singular points. To see the topological property, 
one may pick off these two Dirac points for given param- 
eters. Then, is reduced to — {p2,P2}, which can 
contract to an 5*^ while at the phase boundary it is 
which can contract to origin. The later fact means there 
is a discontinuity of the mapping at the phase boundary. 
Inside of the gapless phase, the Hopf mapping now is re- 
duced to a mapping (u, v) from 5^ to . Taking (px^Py) 
as the coordinate in the source , then the winding num- 
ber of this mapping ^ISi] is the same as that defined in 
two-dimensional Kitaev-type model 

z/ = -3- / dpxdpy{dp^m{p))x{dp m{p))-m{p){7) 

Since {u{oo), v{oo)) — (1, 0) is the north pole of the target 
52 and (w(0),w(0)) = (0,1) is the south pole, we know 
that I' — 1. Therefore, the ground state of the gapless 
phase is topologically non-trivial [l^. ( As explained in 
[ig | , although there is a phase transition from gapless to 
gapped phase transition in Kitaev's original model, both 
phases are topologically trivial.) 

Loops, Knots and Links: We now define the loops, knots 
and links. Loops are excited if Sp — —1. Two loops 
may merge into a single loop by operation of plaquette 
operator. Knot and link excitations are created by this 
operations. To see these operators more clearly, we show 
it by using a model on a multi-layer honeycomb lattice 
which is in the same university class with the diamond 
lattice model in the long wave length limit, i.e., they 
both are described by a 3-d p-wave paired state. The 
links arrangement of the new model is shown in Fig. 2a 
and the Hamiltonian reads 

i 

y\ 1 1+62 1— 1 — e^+ei/ 

+ JyX'^i-es'^i + '^i-e,-e3'^i-eJ> 
+ -'"^;(0-f-e,0-f + f^f-e3^f-e3+eJ> (8) 



FIG. 3: (a) The topological equivalent manifold of model H^, 
a solid torus, (b) A typical link excitation (Hopf link), (c) A 
knot excitation (trefoil). 

where the fundamental translational invariant vector is 
i = m(ei — 63) -|- n(e2 — 63) -f 2Ze2, with basis ei — 

(i,^,0), 62 = (i,-^,0), eg = (-1,0,0), e, - 
(0,0,1). Again, there exist a series of string operators 

— rimJo *"''^mei3 which commutc with the Hamilto- 
nian under periodic boundary condition. Using Majo- 
rana fermions, we obtain bilinear fermionic Hamiltonian 
and the closed string excitation-free Hamiltonian has a 
dispersion E — |£'i£'2|^/^[23]. The gapless phase falls in 
the area given by the inequalities | Jt, | < | | -f | Jy | -|- I Jz | , 
2|J,Jj| < \ Jv? ^(\Jx\ + \Jy\ + \'h\f {i,j=x,y,z). One 
may prove that is also equivalent to a p-wave paired 
model. 

The lattice for is topologically equivalent to a solid 
torus X with e+3 as its tangent vector (Fig. 3a). All 
the loop operators following the geodesic circles commute 
with the Hamiltonian and they are good quantum num- 
bers. Two neighboring loops following zigzag geodesic 
circle are the product of plaquette operators P sand- 
wiched between them, Pi^2----Pn~i = SpS^^^ , where 

Pj — '^j'^2j'^3j'^4j'^5j'^6j the plaquette operator de- 
fined on the hexagon. In this way, we merge two loops 
into a single loop. These loops form a loop gas. Links 
and knots can be obtained by twitted boundary condi- 
tions (See Figs. 3b and 3c). These knots and links may 
also be thought Z2 Wilson loops but they are not exact 
eigen excitations of H^. However, the non-commutativity 
between Hz and the knots and links may only happen in 
the bonds connecting two geodesic circles. Therefore, one 
can ignore this non-commutativity in the thermodynamic 
and long wave length limits and takes these Wilson loops 
to be quasi-exact eigen excitations. The above analysis 
may also be done on diamond lattice. However, since 
there is a shift between two adjacent layers, the illustra- 
tions are not so direct. 

Statistics: One can calculate the linking number of a 
closed string configuration. In the gapped phase, the 
Hopf invariant (O counts the linking number of the dual 
links in momentum space, which is equal to the linking 
number on the original co-ordinate space. The fact that 
the Hopf invariant is the abelian Chern-Sinoms implies 
these closed string excitations obey abelian statistics. 

To see the closed string excitations obey non-abelian 
statistics in the gapless phase, we follow Witten's discus- 
sions in dealing with knots and links in SU(N)fc Chern- 
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FIG. 4: (a) The left and right pieces of cut by (b) 
Three states L+,Lo and L_ of the right piece, which are re- 
lated by the matrix B. 

Simons field theory Q. We compact the 3-d lattice to 
in the continuum limit and then the Wilson loops are 
almost exact eigen excitations. We consider two Wil- 
son loops in S^. Using an S"^ cuts into two pieces 
as shown in Fig. 4a, the left piece may contain very 
complicated stuff but the right piece is simple. Then, 
two Wilson loops puncture four points which are corre- 
sponding to creating two spin fields a and destroying two 
a on S^. The low energy limit of BdG equations ^ is 
the 3-|-l-dimensional Dirac equations. Hence, the Hilbert 
space of the massless Majorana fermions restricted on 5^ 
is determined by c = 1/2 conformal field theory. The fu- 
sion rule a ■ a = I + ijj implies the Hilbert space on this 
punctured is two-dimensional. Therefore, three states 
L_i_,Lo s-nd showed in Fig. 4b are linear-dependent. 
By using the monodroniy B matrix , these states may 
relate to each other through 

Lo = BL+, L^=BLo = B^L+, (9) 

where the B matrix is not proportional to identity. Obvi- 
ously, the operation through the B matrix is an exchange 
of these vortices a and then they obey non-abelian statis- 
tics. Any spin-i has an SU(2) gauge symmetry if the 



spin is represented by fermion or boson operators because 
the Hilbert space is enlarged [13, [HI . So is the present 
model. Fixing Uij = ±1 is corresponding to a gauge 
fixing, i.e., reducing the gauge symmetry from SU(2) to 
Z2. After gauge transformation, a carrys a unit SU(2) 
charge. By using the non-ablelian bosonization [22|, the 
massless Majorana fermion theory on is bosonized to 
an SU(2)2/U(1) coset Wess-Zumino-Witten model. The 
k = 2 turns out that the Wilson loops indeed obey non- 
abelian statistics in the sense Witten defined. Due to 
the gauge group is in its defining representation, the B 
matrix has two eigenvalues Ai = e~*'^/*, A2 ~ —e'-'^^/^. 
Thus, the three states in Fig. 4b satisfy the skein relation 

^ qL+ + _ q-^/^)Lo + q-'l^ - 0, (10) 

with q = e*^/^. Witten took this relation as the definition 
of knot polynomials on 5*^ [i|. As Witten pointed out, 
the physical meaning of the skein relation is clear: all 
links and knots can be resolved by the skein relation. 
Therefore, if we know the physical behavior of unknotted 
and unlinked quantum loop gas, the system with knots 
and links may be understood. This means that the study 
of the quantum loop gas is fundamental fl35 . 

In conclusions, we constructed a 3-d exactly soluble 
spin model with closed string excitations obeying non- 
abelian statistics. This is closely related to Witten's orig- 
inal proposal to the statistics of the Wilson loops in the 
Chern-Simons field theory. A topological phase transi- 
tion between gapless and gapped phases was first time 
predicted. 
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construct another exactly solube model which is on the 
lattice as shown in Fig. 2b. However, the dispersion is 
more complicated. 



